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Abstract 

The dynamical behaviors of two interacting dark energy models are considered. In addition to 
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I. INTRODUCTION 



There exists mounting evidences that the Universe is experiencing accelerated expansion, 
driven by an unknown energy component called "dark energy". The nature and origin 
of dark energy have been an active research topic in the past years. Because the only 
observable effect of dark energy is manifested by gravitational interaction, we know nothing 
about the nature of dark energy except that it has negative pressure. One simple dark 
energy candidate which is consistent with current observations is the cosmological constant. 
However, the small value of the vacuum energy density imposes a big challenge to particle 
physics. Furthermore, the cosmological constant model faces the "coincidence" problem: 
Why is the dark matter energy density comparable to the dark energy density now? 

To alleviate the coincidence problem, other dynamical dark energy models were proposed, 
such as the quintessence model the holographic dark energy model j^, the Chaplygin 
gas model [sl, and the tachyonic model 41. Recently, the weak gravity conjecture was used 
to constrain the property of dark energy 5|. It is also possible that the Einstein theory of 
gravity needs to be modified in order to explain the accelerated expansion. These models 
include the 1/R gravity Q], the f{R) gravity ^j, the DGP model 8], and string inspired 
models [9|]. 

The attractor solution is independent of initial conditions. If the dark energy 
model has an accelerated scaling attractor solution and the ratio of the energy density 
^dark energy /^dark matter betwccn the two dark scctors is Order 1, then the coincidence problem 
can be alleviated. It is well known that for the quintessence model, exponential potentials 
have scaling attractor solutions [l^. For more general scalar field model, scaling attractor 
solutions were also obtained In this Letter, we discuss the dynamical behaviors of the 
quintessence model with exponential potential V^(0) = Voexp(— kA0), here = SnG. Since 
the late-time accelerated scaling attractors of the exponential potential is the scalar field 
dominated solution, fidark energy = 1 [3, it does not provide a satisfactory solution for the 
coincidence problem. In general, the dark energy may not evolve independently, a coupling 
between the dark matter and dark energy is possible 12j. With the interaction between 
the dark sectors, accelerated scaling attractors with Qdark energy /^dark matter = 0(1) can be 
achieved, therefore the coincidence problem can be solved. The interaction between dark 
sectors changes the perturbation dynamics and modifies the cosmic microwave background 
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spectrum 



13|. For the dark energy model with constant equation of motion parameter w, 



it was found that the curvature perturbation has a super-Hubble instabihty in the early 
radiation dominated era, whenever a particular interaction term is present [l^. A more 
careful analysis finds that the stability of the curvature perturbation depends on the form 
of the interaction between dark sectors |l5|. The dynamical quintessence model considered 
in this Letter may not suffer the instability problem. 

By introducing the interaction between dark matter and dark energy, the conservation 
equations become 

Pm + SHpm = -Q, (1) 

+ 3H{1 + w^)p^ = Q, (2) 

where the dark matter energy density is pm, the dark energy density = 0^/2 + V{(f)), the 
dark energy pressure = 0^/2 — V{(t>), the equation of state of the dark energy = p^/prp 
and Q stands for the interaction term. The phenomenological interaction term Q is inspired 
from the interaction between the dilaton field a and the matter field in the scalar-tensor 
theory of gravity 



(3) 



''or a general coupling function ^(cr), the interaction term Q = —3pmH[d(\n^)/d(\na)]/2 



li. 



The Letter is organized as follows. In section 2, we review the method of the phase-plane 
analysis by studying the model discussed in pL8|. In section 3, we discuss the interaction 
model Q = aon'^H'^pf^ and its accelerated scaling attractors. In section 4, we study the 
interaction model Q = j3hi'^H~^pm4>'^ and its accelerated scaling attractors. We conclude the 
Letter in section 5. 



II. INTERACTING MODEL 1 

In this section, we consider the interaction Q = aHp^ [18"] to show the phase-plane 
analysis. Using the dimensionless variables 
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Eqs. ([T]), ([2]) and the Friedmann equation become 

x' = -3x + y^Xy' + h{l + x'- y') + (5) 

y = — ^Ax?/ + -?/(l + a; ), (6) 

where a prime denotes d/dlna. Setting x' = and ?/' = 0, we find that the fixed points of 
the autonomous system ([5]) and are 



{Xci = ±1, Vci = 0), {Xc2 = ±y^' Vci = 0), (xcs = Z/cs = Y 1 - y ), 

a + 3 v^(a + 3)2 - 2aA^ 

^"^^ = TIT' = 713^ ^ 

In terms of the variables x and y, the dark energy density and the equation of state of 
the total matter are 

Q^ = x'^ + y'^, wtot = — V — = a;^ - y'^- (8) 

P0 + Pm 



Since the model was already discussed in [18| and the other fixed points are not interesting for 
our purpose, here we use the fixed point {x^, yci) as an example to discuss the stability of the 
fixed point. For the existence of the the fixed point (xc4, ?/c4), we require (a + S)^ — 2q;A2 > 
and < f2<A < 1. Therefore, we get the existence conditions 



(a + 3) 
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> a + 3, > a > -3 
> A^ > a + 3, 3 > a > 0, (9) 



2a 

To discuss the stability of the fixed point, we need to expand the system ([5]) and ([6]) around 
the fixed point. In general, for an autonomous system 

x' = f{x,y), y' = g{x,y), (10) 

we have a constant nonsingular matrix at the fixed point {xc, yc), 

j^^ = Vc) ai2 = f (Xe, yc) I ^^^^ 

The eigenvalues of the matrix M are 



0-21 = ^{xc, yc) = ^{xc, yc) 



an + a22 ± a/ (an + 022)^ - 4(aiia22 - 012021' 



(12) 



If the real parts of the eigenvalues of the matrix M are negative, then the fixed point is a 
stable point. So the conditions for the fixed point to be stable are 

On + < 0, aiia22 - 012021 > 0. (13) 

Combining equations (j5]), ([6]) and the conditions (fT3l) . we find that the stability conditions 
for the fixed point (xc4, yd) are 

3(a + 3)(a- 1) 



2a 

(a + 3)2 



> > a + 3, > a > -3 

> A^ > a + 3, 3 > a > 0. (14) 



2a 

The result is plotted in Fig. [H From Fig. [H we see that the parameter space for the fixed 
point to be stable is much larger than that obtained in [isl . To verify the correctness of our 
result, we numerically solve the system equations ([5]) and ([6]) with different initial conditions 
for the parameters (a, A) = (0.6, 2.5) and (a, A) = (-0.3, 2.6). The results are shown in Figs. 
[2] and [31 The parameters are outside the parameter space for the fixed point (xc4, ycA) to be 



stable in 



181], but satisfy our conditions ([9]). From Figs. [H [2]and[3l we see that the fixed 



point {xci, yd) is a stable point when (a, A) = (0.6, 2.5) or (a, A) = (-0.3, 2.6). 



III. INTERACTING MODEL 2 

The interaction Q = aHp^ can take the more general form Q = aQK,'^'^~'^ H^""^^ p"^, so 
that the dynamical system still has the attractor solution x"^ + y"^ = 1. In this section, we 
consider n = 2 for simplicity. The autonomous system is 

x' = -3x + + + - y^) + -a^- (15) 



y' = — ^Ax2/ + -?/(l + a;2-2/2). (16) 



For this system, the fixed points are 



(Xcl = ±1, Vcl = 0), (Xc2 = ± J— yc2 = 0), {Xc3 = ?/c3 = A/l " 

1 + ao V6 V 6 



(Xc4, ?/c4 = Y 1 + 2:^4 - y6AXc4/3), (Xc5, l/c5 = Y 1 + 2:^5 - V6AXc5/3), 

where 

_ («o-l)A + VA2(ao-l)' + 12«o 
5 



-5 





X 
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FIG. 1: The region I is the parameter space for the fixed point (xc4, Uca) in model 1 to be stable, 
and 

The accelerated attractors that solve the coincidence problem are the fixed points (xc4, 
l/c4) and {xc5, Ucb)- For the existence of the fixed point (a;c4, yd)) we require 



r- , /3(l + 2ao)2 
y 2q;o(1 + q;o) 

A > v^, -1 < ao < 0, (20) 
A > = A; or - \/3 < A < -A/, ao < -1. 

The stability conditions are 

VS < X < Xu, Oio > 0, 
\/3<A<A+, -l/2<ao<0, 
V3<A<A_, -1< ao < -1/2, 

A; < A < A_, -3 < ao < -1, (21) 




FIG. 2: Phase-space trajectories for the model 1 with a = 0.6 and A = 2.5, the stable fixed point 
is (xc4, yc4)=(0.59, 0.38). 



where 



and 



A, 



-3(1 - 2«o + - 4«o - 28ag - 32«3 
4ao(l + tto) 



3(1 - 2ao - V 1 - 4ao - SSag _ 32al 



(22) 



(23) 



4ao(l + ao) 

The regions of the parameters ao and A for the fixed point (xc4, ?/c4) to be stale are plotted 
in the region I in Fig. HI 

Note that the stability condition for the fixed point (xc3, i/cs) is A^ < 3. From Fig. HI we 
see that for ao < —1 and A < -\/3, both the fixed points (xc3, ycs) and (xc4, yd) are stable 
points, so they are local stable points for those parameters. In other words, when ao < — 1 
and A < -\/3, different initial conditions may lead to either the fixed point {xcs, Vcs) or (xc4, 

To get acceleration, we require Wtot = — y'^ < —1/3. For the stable point (a;c4, yci)-. 
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FIG. 3: Phase-space trajectories for the model 1 with a = —0.3 and A = 2.6, the stable fixed point 
is {xc4, yc4)=(0.42, 0.53). 

the acceleration conditions are 

A > A/, a< -2. (24) 

The accelerated region is also shown in Fig. HI 

For the fixed point (xcs, i/cb), the existence conditions and the stability conditions can be 
obtained from those of the fixed point (xc4, yd) by replacing A with —A in equations (|20|) 
and (12T|) . The existence conditions are 

-V3>X> -Xu, ao>0 

A < -VS, -l<ao< 0, (25) 
> A > A; or A < -A,, ao < -1. 
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FIG. 4: The stability conditions for the fixed points in model 2. The region I is the parameter 
space for the fixed point (xc4, yc4) to be stable and the region II is the parameter space for the 
fixed point {xc5, Vcb) to be stable. The dash-dot line denotes the acceleration condition. 

The stability conditions are 

- \/3 > A > -A„, tto > 0, 
-V3>A>-A+, -l/2<«o<0, 
-V3>A>-A_, -l<«o<-l/2, 
-A/ > A > -A_, -3 < ao < -1, (26) 

The condition (!26l) is shown in the region II in Fig. |H The acceleration conditions are 



A < -Az, a < -2. (27) 

These results are summarized in Table [11 

From Fig. HI we see that the attractors (xc4, Uca) and {xc^, i/cs) lead to accelerated scaling 
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X 


y 


Stability Condition 




Acceleration Condition 


1 





Unstable 


1 


No 


-1 





Unstable 


1 


No 


/ Qo 
Y l+ctQ 





ao > 0, A > A„ 


l+QO 


No 


1 oto 
V l+cto 





oq > 0, A < — Au 


ao 


No 


A/^/6 




A2 <3 


1 


A2 <2 






Equation ^ 




Equation ^ 


Xc5 


yc5 


Equation ([26]) 


^c5 + ^cS 


Equation ^ 



TABLE I: The behaviors of the fixed points in model 2. 



attractors only when -A_ < A < -(4ao/(l + 2ao))^/^ or A_ > A > (4ao/(l + 2ao))^/^ if 
-2 < ao < -1 and -A^ > A > -A_ or A; < A < A_ if -3 < c^o < -2. 



IV. INTERACTING MODEL 3 

In this section, we take the interaction term Q = fSn^^H^''^^ p^cfP' . The dynamical system 
has the attractors x^+y^ = 1. For simplicity, we consider n = 1 case, the autonomous system 
is 

x' = -3x + + + - y^) + 3/3x(l - x^ - y^), (28) 

3 

y' = -—Xxy + -y{l + x^ - y^). (29) 

The fixed points are 



A / A^ 

{xci = ±1, yd = 0), {xc2 = 0, yc2 = 0), (xcs = yc3 = \/ 1 - -^), 

V6 V 6 (30) 

(Xc4, ?/c4 = \ll+xl^- V6ASc4/3), (XcS, l/c5 = ^1 + X^g - V6AXc5/3), 

where 
and 



2v^/3 • ^^^^ 
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For the fixed point (xc3, i/cs), the existence condition is < 6, and the stabihty conditions 

are 

-V6<\<V6, (3> 1/2, 

The condition (133!) is shown in the region I in Fig. [51 The acceleration condition is —\/2 < 
X<V2. 

For the fixed point {xc4, Uca), the existence condition is 
and the stability condition is 

A<-y^, /3<l/2. (35) 

the condition (|35|) is shown in the region II in Fig. [51 The acceleration condition is 

-^3/(1 -p) > A > -y/^ and /? < 0. 

For the fixed point (xc5, Ucb), the existence condition is 

and the stability conditions are 

A>y^, /5<l/2. (37) 

The condition (1371) is shown in the region III in Fig. [51 The acceleration condition is 
v/3/(l - (3) < A < and /? < 0. These results are summarized in Table [Tll 

V. DISCUSSIONS 

We considered two phenomenological interacting models Q = a;o-f^^~^"Pm Q = 
(3H^^^"p'^(j)'^. The scaling attractor solutions x^+y^ = 1 for the non-interacting quintessence 
model with exponential potential remain to be the scaling attractors in the interacting mod- 
els. We have studied the dynamical behaviors of the two interacting dark energy models. 
For the interacting model Q = aoH~^p'^, we find that the fixed points (xcs, i/cs) and (xc4, 
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FIG. 5: The stability conditions for the fixed points (xcs, Ucs), (xca^, yd) and (xcs, ycs) in niodel 3. 
The region I is the parameter space for the fixed point {xcs, Ucs) to be stable, the region II is the 
parameter space for the fixed point {xc4, Uca) to be stable and the region III is the parameter space 
for the fixed point {xc^, ycs) to be stable. The dash-dot line denotes the acceleration condition. 

Uci) are stable points if ao < —1 and A < -\/3, and the fixed points (xcs, ya) and (a;c5, Z/cs) 
are stable points if ao < —1 and A > — \/3- In other words, in the parameter region ao < — 1 
and A < or A > — -\/3, the fixed points (xc3, ^cs), {xca-, Uca) and (xc5, Ucb) are local stable 
points. This type of local stable points are new in the dark energy models. 

These models have the late time accelerated scaling attractors with Vt^/Vt^ = 0(1). We 
can easily match Q^/Q^ to observations by a simple choice of parameters. Since the solution 
is a scaling attractor, the value of fi^^/fi^ is insensitive to initial conditions, therefore the 
why now problem is resolved. For the interacting model Q = f3H~^ pmcf^ ■, when we choose 
[3 = —2.0 and A = 2.5, the stable fixed point is (xcs, yc5) = (0.31, 0.68) with wtot = —0.37 
and Qfp = 0.75. Note that to get accelerated attractor solution and alleviate the coincidence 
problem, we find that a < and (3 < 0, so the energy transfer goes from dark energy to 
dark matter. This result is easily understood. The energy transfer from dark energy to 
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X 


y 


Stability Condition 


^0 


Acceleration Condition 


1 





/5 > 1/2 and A > 


1 


No 


-1 





/3 > 1/2 and A < -V6 


1 


No 








Unstable 





No 




x/l - AV6 


Equation ([33|) 


1 


A2 <2 






Equation ([35]) 




A > and /? < 




2/c5 


Equation ([37|) 


^c5 + Vcb 


A < and ^ < 



TABLE II: The behaviors of the fixed points in model 3. 



dark matter makes the dark matter to decrease slower and dark energy to decrease faster, 
therefore alleviating the coincidence problem. 
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